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Abstract 

In this work, we extend the non-orthogonal amplify-and-forward (NAF) cooperative diversity scheme 
to the MIMO channel. A family of space-time block codes for a half-duplex MIMO NAF fading 
cooperative channel with N relays is constructed. The code construction is based on the non-vanishing 
determinant (NVD) criterion and is shown to achieve the optimal diversity-multiplexing tradeoff (DMT) 
of the channel. We provide a general explicit algebraic construction, followed by some examples. In 
particular, in the single-relay case, it is proved that the Golden code and the 4x4 Perfect code are optimal 

(N ■ 

;> for the single-antenna and two-antenna case, respectively. Simulation results reveal that a significant 

gain (up to 10 dB) can be obtained with the proposed codes, especially in the single-antenna case. 
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Cooperative diversity, relay channel, amplify-and-forward (AF), multiple-input multiple-output (MIMO), 
space-time block code, non-vanishing determinant (NVD), diversity-multiplexing tradeoff (DMT). 
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H NOTATIONS 



In this paper, we use boldface lower case letters v to denote vectors, boldface capital letters 
M to denote matrices. CM represents the complex Gaussian random variable. E[-] stands for the 
expectation operation and [-] T , [•]* denote the matrix transposition and conjugated transposition 
operations. ||-|| is the Euclidean vector norm and ||-|| F is the Frobenius matrix norm. |«S| is the 
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cardinality of the set S. (x) + means max(0, x). R, C, Q and Z stand for the real field, complex 
field, rational field and the integer ring respectively. For any quantity q, 

q = SNR a means lim - lo f^ = a 

SNR^oologSNR 

and similarly for < and > . 

I. Introduction 

On a wireless channel, diversity techniques are used to combat channel fadings. Recently, there 
has been a growing interest in the so called cooperative diversity techniques, where multiple 
terminals in a network cooperate to form a virtual antenna array in order to exploit spatial 
diversity in a distributed fashion. In this manner, spatial diversity gain can be obtained even 
when a local antenna array is not available. Since the work of [1], [2], several cooperative 
transmission protocols have been proposed [3]-[9]. These protocols can be categorized into 
two principal classes : the amplify-and-forward (AF) scheme and the decode-and-forward (DF) 
scheme. In practice, the AF scheme is more attractive for its low complexity since the cooperative 
terminals (relays) simply forward the signal and do not decode it. Actually, for most ad hoc 
wireless networks, it is not realistic for other terminals to decode the signal from a certain user, 
because the codebook is seldom available and the decoding complexity is unacceptable in most 
cases. 

The non-orthogonal amplify-and-forward (NAF) scheme was proposed by Nabar et al. [5] for 
the single-relay channel and was then generalized to the multiple-relay case by Azarian et al. 
[6]. In [6], it is shown that the NAF scheme outperforms all previously proposed AF schemes in 
terms of the fundamental diversity-multiplexing tradeoff (DMT) [10] and that it is optimal within 
the class of AF schemes in the single-relay case. The superiority of the NAF scheme comes 
from the fact that the source terminal is allowed to transmit during all the time, which boosts up 
the multiplexing gain. However, even though they showed that the DMT of this scheme can be 
achieved using a Gaussian random code of sufficiently large block length, no practical coding 
scheme that achieves the tradeoff has been proposed since then. 

The main contributions of our work are summarized as follows: 

1) We extend the single-antenna NAF scheme proposed in [5], [6] to the multiple-antenna case. 
We establish a lower bound on the optimal DMT of the MIMO NAF channel. In particular, 
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we show that the maximum diversity order of a single-relay MEMO NAF channel is lower- 
bounded by the sum of the maximum diversity order of the source-destination channel and 
the maximum diversity orders of the source-relay-destination product channels. This lower 
bound is tight when the source, relay and destination antenna number n s , n r and n d satisfy 



2) We provide an explicit algebraic construction of short block codes that achieves the optimal 
DMT of the general multiple- antenna multiple-relay NAF cooperative diversity scheme. Our 
algebraic code construction is inspired by the non-vanishing determinant (NVD) space-time 
codes design for MIMO Rayleigh channels [11]. First, we show that for any linear fading 
Gaussian channel (not only the Rayleigh channel as in [10], [11]) 



in the high SNR regime, the error event of a "good" space-time code X (which will be 
properly defined later) occurs only when the channel is in outage. Therefore, the optimal 
DMT can always be achieved by X. Then, as in [6], we derive equivalent signal models 
of the AF cooperative schemes in the form (HI, subject to certain input constraint (e.g., 
block diagonal for multiple-relay channel). Since codes that achieve the optimal DMT of 
the equivalent channel (HJ) also achieve the optimal DMT of the corresponding cooperative 
channel, optimal codes for an AF cooperative channel can be obtained from the NVD 
criterion. As a result, we show that for a single-relay AF channel with n s antennas at the 
source terminal, a 2n s x 2n s full rate NVD space-time code (e.g., the Golden code [12] for 
the single- antenna case and the 4x4 Perfect code [13] for the two-antenna case) can be 
directly applied to construct an optimal block code. In the iV-relay case, the optimal code is 
constructed from a block-diagonal NVD space-time code with N blocks. The performance 
of our construction is confirmed by simulation results. 
The rest of the paper is outlined as follows. Section |n| introduces the system model and 
recalls the single- antenna NAF protocol as well as the equivalent channel model. In section [Oil 
we extend the NAF scheme to the MIMO channel and develop a lower bound on the optimal 
diversity-multiplexing tradeoff. The codes design criteria are derived in section|lV|and the explicit 
algebraic construction that satisfies the design criteria is provided in section |Vj Section IVTI shows 



n s — n&\ >7h 



y = VSNRHx + z, 



(1) 
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some examples of channel configuration and the parameters of the corresponding optimal codes. 
Simulation results on our construction are available in section IVIII Section IVIIII contains some 
concluding remarks. For continuity of demonstration, most proofs are left in the appendices. 

II. System Model and Problem Formulation 

A. Channel Model 

We consider a wireless network with N + 1 sources (users) and only one destination. The 
channels are slow fading (or delay-limited), i.e., the channel coherence time is much larger than 
the maximum delay that can be tolerated by the application. For the moment, we assume that 
all the terminals are equipped with only one antenna. The multi-antenna case will be treated 
separately in section |ni| The channel is shared in a TDMA manner, i.e., each user is allocated a 
time slot for the transmission of its own data. Within the same time slot, any of the other N users 
can help the current user transmit its information. The extension to a more general orthogonal 
access scheme is straightforward. Suppose that the network configuration is symmetric. Without 
loss of generality, we consider only one time slot and the channel model becomes a single-user 
relay channel with one source, N relays and one destination, as shown in Fig. [U Here, we 
exclude the multi-user case, where information of more than one user can circulate at the same 
time in the network {e.g., the CMA-NAF scheme proposed in [6]). 

s / d 




o 



Fig. 1. A relay channel with one source (s), one destination (d) and TV relays (Ti, . . . , fjv). 

In Fig.[H variables /, hi and g^i = 1, . . . , N stand for the channel coefficients that remain 
constant during a block of length L. As in the previous works we cite here, we assume that 
all the terminals work in half duplex mode, i.e., they cannot receive and transmit at the same 
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time. The channel state information (CSI) is supposed to be known to the receiver but not to 
the transmitter. 

B. The Non-Orthogonal Amplify '-and-Forward Relay Channel 

In our work, we consider the NAF protocol ( [5], [6]). In this scheme, the relays simply scale 
and forward the received signal. However, unlike the orthogonal AF protocols, the source can 
keep transmitting during the transmission of the relays. 





x\ 




by] 



I I transmitted 
L _ _ _ _~J received 



Fig. 2. The NAF frame structure of a single-relay channel, b is the normalization factor such that by r is subject to the power 
constraint. 

1 ) The Single-Relay Case: In the single-relay case, each frame is composed of two partitions 
of T/2 symbols 1 . The frame length T is supposed to be smaller than the channel coherence 
time L, i.e., the channel is static during the transmission of a frame. The half duplex constraint 
imposes that the relay can only transmit in the second partition. The frame structure is illustrated 
in Fig. El from which we get the following signal model 

Vv = VtmSNR/x, + v v 

< y r = a/ttipSNR/ixi + w (2) 

k y 2 = y/n 3 SNRg (by r ) + y^SHRfx 2 + v 2 

where x,. y, 6 C T//2 , % = 1, 2 are the transmitted signals from the source with normalized power 
and the received signals at the destination, respectively, in the i th partition; y r 6 C T ^ 2 is the 
received signal at the relay in the first partition; V\.v>.w E C T//2 are independent additive white 
Gaussian noise (AWGN) vectors with i.i.d. unit variance entries; the channels between different 
nodes are independently Rayleigh distributed, i.e., f,g,h ~ £/V(0, 1); p is the geometric gain 

'it is shown in [6] that giving the same length to the two partitions is optimal in terms of the DMT. 
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Fig. 3. The NAF frame structure of an iV-relay channel. Solid box for transmitted signal and dashed box for received signal. 



representing the ratio between the path loss of the source-relay link and the source-destination 
link; b is the normalization factor satisfying E {||6j/ r || 2 } < \, i.e., 

|6| 2 < 1 — 2 . 

~ 7ripSNR|fc| 2 + 1 

We consider a short term power constraint, i.e., the power allocation factors 7Tj's do not depend 
on the instantaneous channel realization f,g and h, but can depend on p and SNR. We impose 
that J2i n i = 2 so that SNR denotes the average received SNR at the destination 2 . 

As shown in [6], the channel model © is equivalent to T/2 channel uses of a 2 x 2 channel : 



V^SNR/ 







Xi + 







W; 



Vi for i = T/2 



where «j = [«i[i] w 2 [i]] T for m e {x,y,t;} and Uk[i] denotes the i th symbol in the k th partition. 
In the following, we consider a more convenient normalized model 

Vi = VSNRHxi + Zi for i = 1, . . . , T/2 

where Zi ~ CAf(0, 1) is the equivalent AWGN and 

/ 5rT/_ 



H 



7T1 7T3 pS N R 

l+7r 3 SNR|6 5 | : 



l+7r 3 SNR|6 ff |" 



(3) 



2 The total transmit power in two partitions are (7Ti + 7T2 + 7T3)SNR. Since the channel coefficients and the AWGN are 
normalized, (7Ti + 7T2 + 7T3)SNR represents the average received SNR per two partitions as well. 
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2) The Multiple-Relay Case: In the multiple-relay case, a superframe of iV consecutive 
cooperation frames is defined, as shown in Fig.|3] It is assumed that the channel is static during 
the transmission of the whole superframe (of NT symbols). The N relays take turns to cooperate 
with the source. Within each cooperation frame, the cooperation is in exactly the same manner 
as in the single-relay case. However, by allowing an encoding over the whole superframe, a 
diversity order of N + 1 is achieved. 

C. Diversity-Multiplexing tradeoff (DMT) 

Definition 1 (Multiplexing and diversity gain [10]): A coding scheme {C(SNR)} is said to 

achieve multiplexing gain r and diversity gain d if 

fl(SNR) logP e (SNR) 
hm - — — — = r and hm — - — — — — = —a 

SNR^oo logSNR SNR->oc logSNR 

where i2(SNR) is the data rate measured by bits per channel use (PCU) and P e (SNR) is the 
average error probability using the maximum likelihood (ML) decoder. 
The optimal DMT of the single-antenna iV-relay NAF channel 

cW(r) = (1 - r)+ + N(l - 2r)+ (4) 

is found in [6], where the achievability is proved by using a Gaussian random code with a 
sufficiently long block length. 

III. The NAF Scheme for MIMO Channel 

In this section, we generalize the NAF protocol to the MIMO case, where each terminal is 
equipped with multiple antennas. The notation (n s ,n n n d ) will be used to denote a single-relay 
channel with n s , n r and antennas at the source, relay and destination. All matrix variables 
defined for a single-relay channel apply to a multiple-relay channel with an index i denoting the 
i 1 * 1 relay. 

A. Signal Model 

For convenience of demonstration, we only present the signal model of the single-relay 
channel. The extension to the multiple-relay case is straightforward. 
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1) n r < n s : In the case n T < n s , a direct generalization of © is as follows : 



y, = y / n 1 SNRFX 1 +V 1 

Y, = ^/ti iP SNRHX 1 + W 

Y 2 = y/n 3 SmG(BY r ) + y/7r 2 SMRFX 2 + V 2 



(5) 



where F, G, H are n A x n s , n d x n r and n r x n s independent matrices, respectively, with zero mean 
unit variance i.i.d. Gaussian entries; Xj's are n s x -| matrices with i.i.d. zero mean unit variance 
entries, representing the space-time signal from the source; V \.V> and W are independent 
AWGN matrices with normalized i.i.d. entries; the power allocation factors 7T;'s satisfy n s {jii + 
7r 2 ) + n r ii 3 = 2 so that SNR denotes the received SNR per receive antenna at the destination; B 
is an n r x n r matrix equivalent to the "normalization factor" b in the single- antenna case and is 
subject to the power constraint E| H-BYVHp j < ^n r which can be simplified to 

Tr|(I + mpSNRHW) BtB} < n r . (6) 

Now, as in the single-antenna case, we obtain an equivalent single-user MIMO channel 



y, = VSNRHxi + z t for % = 1, . . . , T/2 (7) 

where X{ = \X\ [i] T X > [i] T ] T and jjfj = fYi [i} T Y 2 [i] T ] T are the vectorized transmitted and received 
signals with M[i) denoting the i th column of the matrix M; z,, ~ CjV(0, I) is the equivalent 
AWGN; the equivalent channel matrix H is 

y/^lF 



H 



(8) 



with 



P = VSNRGB 



(9) 



and T being the whitening matrix satisfying (r f r) _1 = (IT 1 ) -1 = I + n 3 PP' ! = E. 

2) n r > n s : With more antennas at the relay than at the source, the relay can do better than 
simple forwarding. In this case, the received signal at the relay is in the n s -dimensional subspace 
generated by the n s eigenmodes of H, represented by Uh from the singular value decomposition 
H = U h^hV^. Since the relay-destination channel G is isotropic, it is of no use to forward 
the received signal in more than n r antennas (spatial directions). 
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However, by using only a subset of the antennas at the relay, we cannot obtain all the diversity 
gain provided by the channel G. To exploit all the available diversity, we propose two schemes. 
The first scheme is the virtual iV-relay scheme. Since no CSI is available at the transmitter 
and therefore no antenna combination is a priori better than the others, one solution is to use 
all the (™ r ) antenna combinations equally. Intuitively, the typical outage event is that all the 
antenna combinations are in deep fade, which implies that the channel G is also in outage. In 
this scheme, a superframe of (™ r ) cooperation frames is constructed. Within each cooperation 
frame, a different combination of n s antennas is used. The second scheme is the antenna selection 
scheme, which is used only when limited feedback from the destination to the relay is available. 
In this case, the destination tells the relay which antenna combination is optimal according 
to a given criterion (e.g., maximization of mutual information). With this scheme, maximum 
diversity gain is obtained without a superframe structure, which means a significant reduction 
of coding-decoding complexity. 

As an example, let us consider a (1, n T , 1) relay channel with n T > 1. In this case, the received 
signal at the relay can be projected into a one-dimensional subspace. Consider a superframe of 
n r cooperation frames. In each cooperation frame, a different relay antenna is used. This scheme 
is virtually an iV-relay single- antenna channel. The only difference is that the equivalent source- 
relay link, which is \J~^2~\h^ after the matched filter operation, is the same for the N virtual 
relays. In this scheme, the achievable diversity order of the source-relay-destination link is n r , 
since the channel is in outage only when \J^~\h^ is in deep fade or all the n r relay-destination 
links are in deep fade. When feedback is possible, the antenna selection scheme can be used. 
The difference from the first scheme is that only the antenna with maximum relay-destination 
channel gain (say, |<? max | — max is used in the relaying phase. Therefore, only one 

i=l. ..n r 

cooperation frame is needed and the diversity order is also n r . 

B. Optimal Diversity-Multiplexing Tradeoff: a Lower Bound 

With the discussion above, considering the case n r < n s is without loss of generality. In 
addition, since the destination is usually equipped with more antennas than the relays are in 
practice, we will restrict ourselves to the case n d > n r hereafter. In the rest of this section, we 
study the optimal DMT of the MIMO NAF cooperative channel. Unlike the single- antenna case, 
a closed form expression of the DMT of the MIMO NAF channel is difficult to obtain, since the 
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probability distribution (in the high SNR regime or not) of the eigenvalues of H defined in ® 
is unknown. In the following, we will derive a lower bound on the tradeoff, as a generalization 
of the DMT of the single- antenna NAF channel provided in [6]. To this end, we first study the 
DMT of a Rayleigh product channel. 
1 ) DMT of a Rayleigh Product Channel: 

Proposition 1: Let G, H be n x 1,1 x m independent matrices with i.i.d. entries distributed as 
CJ\f(Q, 1). Assume that m > l,n > I and define A = GH, A = \m — n\ , q = min {m, n}, then 
the optimal DMT curve d* A {r) of the Rayleigh product channel, i.e., the channel defined by 



/SNR , 

y = \ -, Ax I z (10) 

V t ■ m 

with z ~ CjV(0, 1) being the AWGN is a piecewise-linear function connecting the points 

(s, d A (s)), s = 0, . . . , I, where 



dl(s) = (l-s)(q-s)-^ 



[(I -A 



2 



(ID 

where a Rayleigh product 



2 

Proof: See AppendixlTTI A more general result is given by [20], 
channel is seen as a special case of the double scattering channels and the assumption m > 
l,n > I is unnecessary. ■ 
Remark 1: From propositionHJ we note that 

(i) d* A (r) only depends on A and q, which means that interchanging m and n does not change 
d* A (r), which is obvious if we consider the fact that AA f has the same eigenvalues as A ] A; 

(ii) d* A (s) is upper-bounded by cIa(s) = min{ dn (s), dc(s)} = (I — s)(q — s) and coincides 
with it when s > I — A — 1; 

(iii) when A > I — 1, d A (s) = Ja(s),Vs. Without loss of generality, assume that n > m. 
Intuitively, by increasing n and keeping m unchanged, the diversity gain of the channel 
G is increasing and so is that of the product channel A. When n — m = I — 1, the best 
possible DMT for q = m is achieved. With this n, the fading effect of G vanishes, since 
d A (s) = dn(s). In other words, it is of no use to have n > I + m — 1, for m and / fixed. 

The optimal DMTs of a Rayleigh product channel for / = 2,3 are illustrated in Fig.|H As 
indicated in remark 1, when A* = 1, 2 for I = 2,3, the tradeoffs of the Rayleigh product channels 
are the same as those of the corresponding Rayleigh channels, i.e., the 2x2 and 3x3 Rayleigh 
channels. 
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2) DMT of a MIMO NAF Channel- 
Theorem 1: For a single-relay MIMO NAF channel ©, we have 

cW(r) > d F (r) + d GH (2r) (12) 

where dp(r) and doH(r) are the optimal DMT of the fading channel F and the fading product 
channel GH, respectively. 

Proof: See AppendixUnl ■ 

The interpretation of (fT2l is as follows. First, since the transmitted signal passes through the 
source-destination link all the time, a diversity gain dp(r) can be obtained. Then, due to the half 
duplex constraint, only half of the transmitted signal is protected by the source-relay-destination 
link, i.e., the channel defined by GH. Fig. |5] shows the lower bound (fT2l for a Rayleigh channel 
with n s = n r = n<j = 2 and 3. 

This result of theorem [T] can be generalized to the iV-relay case. Let Gi and Hi be channels 
related to the i th relay and be similarly defined as G and H . The following theorem gives a 
lower bound on the optimal DMT of the iV-relay MIMO NAF channel. 
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Fig. 5. Lower bound on the optimal DMT: n s = n, = m = 2, 3. 



Theorem 2: For an iV-relay MIMO NAF channel, we have 

N 

d NAP (r) > d F {r) + min ^^(2%). 
In particular, when all the relays have the same number of antennas, we have 

<W0 >d F {r)+Nd GH {2r). (13) 



Proof: See Appendix|IVl ■ 
From ©, we see that the bound in ([131) is actually the optimal tradeoff in the single-antenna 

case {i.e., n s = n r = n d = 1). 

Corollary 1: Let d{ = min{dGi(0), cf^O)}, i = 1, . . . , AT. Then, we have 



TV AT 



<fr(0) ^ ^naf ( ) > d P (0) +X) c Wr i (0). (14) 

i=l i=l 

If all the channels are Rayleigh distributed and \m — n\ > li — 1, Vi, then, we have 



iV 



d NAF (0) = d F (0) + J2d l (0). (15) 

i=l 

Proof: See AppendixQAT] ■ 
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IV. Optimal Codes Design Criteria 

In this section, we will derive design criteria for a family of short codes to achieve the optimal 
DMT of an iV-relay MIMO NAF channel. 

A. A General Result 

Let us first define the "good" code mentioned in section Q] 

Definition 2 (Rate-n NVD code): Let A be an alphabet that is scalably dense, i.e., for < 

r < n 



Then, an n T x nx space-time code X is called a rate-n NVD code if it 

1) is .4- linear 3 ; 

2) transmits on average n symbols PCU from the signal constellation A; 

3) has the non-vanishing determinant (NVD) property 4 . 
The following theorem is fundamental to our construction. 
Theorem 3: For any linear block fading channel 

y = VSNRHx + z 

where H is an n R x n T channel and z ~ £/V(0, 1) is the AWGN, the achievable DMT of a 
rate-n NVD code X satisfies 



where q = min{nR,n T } and d out (r) is the outage upper bound of the DMT for the channel H. 



In particular, for a full rate code (n = q), the upper bound d out (r) is achievable. This theorem 
implies that the NVD property is fundamental for X to achieve all the diversity gain d, for any 
linear fading channel. For a given diversity gain d, the achievable multiplexing gain r of such 
X is a shrunk version of r out (d), the best that we can have for channel H. 

3 X is ^4-linear means that each entry of any codeword X G X is a linear combination of symbols from A. 
4 NVD means that \det(Xi - Xj)\ > k > 0, V'X,.Xj £ X,Xi ^ Xj with k a constant independent of the SNR. 



|^4(SNR)| = SNR™ and 
a e A(SNR) |a| 2 <SNR- 




(16) 



Proof: See AppendixM 
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One of the consequences of theorem |3] is the possibility of constructing optimal codes (in 
terms of the DMT) based on the NVD criteria for some channels. For example, we can get 
an equivalent MIMO space-time model for the single- antenna fast fading channel (also called a 
Gaussian parallel channel) as 



Y = (liag(/ M h N )X + Z 



(17) 



with X, Y and Z diagonal Nx N matrices. The best code that we can have is a rate-1 NVD code 
X due to the diagonal constraint. According to theorem|3l we have dxif) > d 0Ut (Nr) with d out (r) 
the DMT of (flTt without the diagonal constraint. In fact, we can verify that d out (Nr) coincides 
with the DMT of the fast fading channel. The NVD criterion includes the product distance 
criterion since the determinant of a diagonal matrix is the product of the diagonal entries. In 
addition, it implies that the product distance should be non-vanishing as the constellation size 
increases. 

Note that another such general result as theorem has been derived independently in [14]. In 
[11], the NVD property is derived from the mismatched eigenvalue bound (worst case rotation) 
while the results in [14] are derived using the worst case codeword error probability, which is 
effectively the same thing as the worst case rotation. Theorem |3] is a generalization of the result 
in [11] (for the full rate codes) to a rate-n code. This result is more adapted to the algebraic 
construction of explicit codes for the relay channel. 

B. Design Criteria 

With theorem[3l we are ready to give out the design criteria of the optimal codes for the NAF 
cooperative channel. The following theorem states the main result of our work. 
Theorem 4: Let X be a rate-(2n s ) NVD block diagonal code, i.e., 

■ 



X 



— 1 



VX G x 



■■• E N 

where Ej's are 2n s x 2n s matrices. Now consider an equivalent code C whose codewords are in 
the form 



Ci 



c 



N 
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with 



Ci= Sj (l:n S) l:2n s ) Hj (n s +l:2n s , 1 :2n s ) 



Then, C achieves the optimal DMT of the iV-relay MIMO NAF channel with n s transmit antennas 
at the source, by transmitting Cj in the i th cooperation frame. The code C is of length ANn s . 



In section [Tnl a lower bound on the optimal DMT of a MIMO NAF channel is derived. Here, 
theorem |4] shows that the exact optimal tradeoff can always be achieved by a code C, even 
though we cannot obtain its closed form expression. 



A. Notations and Assumptions 

We assume that the modulation used by the source is either a QAM or a HEX modulation. 
The fields representing the modulated symbols will be either Q(i) or Q(j). We denote it as P. 
For each algebraic number field K, the ring of integers is denoted Ok- 

B. Behavior of the Codewords 

We recall that a codeword X is represented by a block diagonal matrix 



with Hj, i = 1, . . . , N being a square 2n s x 2n s matrix. The criteria to fulfill are the following : 

1) full rate : the number of QAM or HEX independent symbols in a codeword is equal to 
iV • (2n s ) 2 corresponding to a multiplexing gain of n s symbols PCU; 

2) full rank : 



Proof: See AppendixlVIl 



V. A Unified Construction Framework 







X 



(18) 








min rank(Xi - X 2 ) = N ■ (2n,); 



(19) 



3) non-vanishing determinant : 



min |det(Xi — X 2 )\ 2 > k 



(20) 



x 1 ,x 2 ec 



with k being some strictly positive constant. 
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C Codes Construction 



We use the same methods as in [13]. Some particular cases can be found in [12], [15]. The 
main difference is in the choice of the base field F. In [13], this base field was equal to P. Here, 
we choose a Galois extension of P with degree iV and denote Tj, i = 1, . . . , N the elements 
of its Galois group Qal-g/f. Now, we construct a cyclic algebra whose center is F. We need a 
cyclic extension over F of degree 2n s . We denote it K. The generator of its Galois group is a. 
The code construction needs two steps. 

1) Construction of the cyclic algebra 



A 



2n s -l 



zt ■ e* 



i=0 



Zi e K 



(21) 



such that e 2ns = 7 G F and ■ e = e ■ a 



). In the matrix representation, we have 

1 ■••0 

001 •. ; 
; ■-. ■•. 
; ■•. 1 
7 ■•• 

2) Application of the embeddings of F/P. 
In terms of matrices, we construct, in step 1, the square 2n s x 2n s matrix S. Then, by applying 
the embeddings of F/P, the codeword is 





and Zi = diag {z h a (z;) , a 2 



X 



(22) 



••• r N (E) 

where we can identify Tj(H) = from dl8b . As usual, we restrict the information symbols to 
be in C P , that is, Z[i] (QAM symbols) or Z[j] (HEX symbols). So, instead of being in F, we 
will be in Of and in the same way, we will be in Ok instead of K. The infinite space-time code 
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is defined as being the set of all matrices 

'2n s -l 

• i E 



C = < 



X 



i=0 




(23) 



D. Codes Properties 

Lemma 1: The code C of (E3l is full rate. 

Proof: In the submatrix H, there are 2n s independent elements of Ok- Each element in 
Ok is a linear combination of 2n s elements of ¥ . Finally, each element of ¥ is a linear 
combination of iV QAM or HEX symbols. So, each codeword X is a linear combination of 
N ■ (2n s ) 2 QAM or HEX symbols. ■ 



Lemma 2: If 7, 7 2 



-7 



2n s -l 



^ Nk/f(K), then the code C is full rank. 



Proof: In [16], it is proved that if 7,7 2 



■ 7 



2n 3 -l 



^ A/k/f(K), then the cyclic algebra A 



is a division algebra (each element has an inverse). 



Lemma 3: If 7, 7 2 , 
more precisely 



.7 



2n,-l 



^ 7Vk/f(K), then the code C has a non-vanishing determinant, 



= min|detXr e Z+\{0} > 1. 
Proof: Because of the structure of X, its determinant is 



N 



N 



detX = H det 7i(E) = J] r, (det(S)) . 



t=i 



But, det(S) is the reduced norm of 2^=0 1 z i e% mus ^ belongs to Of. So, 



A? 



JJri (det(S)) = N ¥/¥ (det(S)) G 



i=i 



with C P = Z[i] or C P = Z[?']. Since det(X) ^ unless X = 0, we get <5 min > 1. ■ 
Finally, the following result is derived. 

Theorem 5: The code C of (l23l) with £j e Ok or a subspace of Ok (which will be in the 
following an ideal of Ok) achieves the DMT of the MIMO NAF cooperative channel when iV 
is the number of relays and n s is the number of antennas at the source. 

Proof: The proof is straightforward and uses the results of the 3 above lemmas. ■ 
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E, Shaping 

As in [12], [13], we may be interested in constructing codes that achieve the DMT and that 
behave well in terms of error probability even for small alphabets such as QPSK (4QAM). In 
that case, we add another constraint to our codes design, the shaping factor. This new constraint 
implies that |7| = 1. Moreover, as in [12], [13], the linear transform that sends the vector 
composed by the iV • (2n s ) 2 QAM or HEX information symbols to vec(X) has to be unitary. 
The following examples will illustrate this claim. 

VI. Some Examples 

We give some examples of the code construction. Our code for an iV-relay A;-antenna channel 
is denoted C Ntk . 

A. The Golden Code [12] is Optimal for the Single-Relay Single-Antenna NAF Channel 

In the case of single-relay single- antenna channel, the codewords are 2x2 matrices. Because 
the Golden code satisfies to all the criteria of subsection IV-B1 it achieves the optimal DMT of 
the channel. 

B. Two Relays, Single Antenna 

Optimal codes for the case N > 1 relays cannot be found in the literature. 

For the 2-relay case, we propose the following code. Codewords are block diagonal matrices 
with 2 blocks. Each block is a 2 x 2 matrix. Let P = Q(z) and F = Q (Cs) with Cs = be 
an extension of Q(i) of degree 2. We choose K = F (Vo) = Q (Cs> In f act > we tr y to 
construct the Golden code on the base field Q (Cs) instead of the base field Q(i). Moreover, the 
number 7 is no more equal to i because i is a norm in Q (£ 8 ) (i = iV K/ / F (Cs))- We choose here, in 
order to preserve the shaping of the code, 7 = Cs- We prove in appendix IVIII that Cs ]Vk/f(K) 
and thus that this code satisfies to the full rank and the NVD conditions. Such a code uses 8 
QAM symbols. Let 6 = ^ ± ^-, 7 = Cs and a : 6 i-> 6 = The ring of integers of K is 

O k = {a + b8 I a, 6 e Z [Cs]}- Let a = 1 + i — iO and a = 1 + i — i6. Codewords are given by 

[h 

X = 

r(H) 
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with 

„ _ 1 a • (si + s 2 Cs + s 3 9 + s 4 C 8 6») a • (s 5 + s 6 ( s + s 7 9 + s 8 ( 8 8) 

VE ( 8 « . ( S5 + S6 ( 8 + S7 # + Ss ( s e) a-(si + s 2 ( 8 + s 3 6 + s 4 ( 8 6) 

and r changes Cs into — Cs- 



C. Four Relays, Single Antenna 

The generalization to N = 4 relays is straightforward. Codewords are block diagonal matrices 
with 4 blocks. Each block is a 2 x 2 matrix. Let P = Q(i) and F = Q with £i6 = be an 
extension of Q(i) of degree 4. We choose K = F (V5) = Q (Ci6, V^)- We choose here, in order 
to preserve the shaping of the code, 7 = £i6- We prove in appendix IVIIII that Ci6 ^ -Wk/f(K) 
and thus that this code satisfies to the full rank and the NVD conditions. 

D. Single Relay, Two Antennas 

Since n s = 2 and N — 1, we need a code whose codewords are represented by a 4 x 4 NVD 
space-time code. The 4x4 Perfect code of [13] satisfies to all criteria. 

E. Two Relays, Two Antennas 

We assume here that the source uses 2 antennas and that there are 2 relays. The idea is to 
construct a 4 x 4 Perfect code not on the base field F = Q(i) as it is the case in [13], but on 
the base field F = Q (£ 8 ). Thus, a rate-2 NVD code can be constructed as follows: 

. TakeF = Q(C 8 ). 

. Choose K = F(2cos(ff)) = Q 2 cos (ff )) . 
• Finally, take 7 = Cs- 

We can show, in the same way as in appendix IVII1 that if ( 8 was a norm in Q (Cs, 2 cos (ff )), 
then i must be a norm in Q (i, 2 cos (ff )) which contradicts the results of [13]. The case of 7 2 
is obvious since 7 2 = i. Now, in order to prove that 7 s = i($ is not a norm, it is enough to 
replace ( 8 by i( 8 and z by y in appendix IVIII and show in the same way that if i( 8 was a norm 
in Q (C 8 ,2cos (ff)), then i must be a norm in Q (i,2cos (ff)). 
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VII. Numerical Results 

In this section, we provide the simulation results on the performance of some of the codes 
proposed in section |VI| The performance is measured by the frame error rate (FER) vs. receive 
SNR per bit. For simplicity, we set the power allocation factors t\\ = 2n 2 = 2tt 3 for all the 
scenarios that we considered in this section. An optimization on the 7Tj's in function of p and 
SNR can improve the performance 5 . However, this kind of optimization is out of the scope of this 
paper and will not be considered here. The transmitted signal constellation is 4- and 64-QAM. 
The geometric gain p varies from to 20 dB. 

A. Single-Antenna Channel 

Fig. [7] shows the performance of the Golden code on the single-relay single-antenna channel. 
The performance of the channel without relay is also shown in the figures. In this case, the frame 
length is 4 symbols. Compared to the non-cooperative case, the Golden code achieves diversity 2. 
For 4-QAM, a gain of 12.5 dB (resp. 13.8, 14.3 and 14.8 dB) is observed for p = OdB (resp. 5, 10 
and 20 dB) at FER = 10 -4 . First of all, note that in the low SNR regime, the non-cooperative 
channel is better than the cooperative channel. This is due to the error cumulation (at the relay) 
which is more significant than the diversity gain provided by the relay in this regime. Then, 
we see that the difference between p = 10 dB and p = 20 dB is negligible, which means that 
a geometric gain of lOdB is enough to achieve the (almost) best performance of the Golden 
code. In practice, it is often possible to find this kind of "helping agent" (with a geometric 
gain of 10 dB). When we increase the spectral efficiency (64-QAM), same phenomena can be 
observed except that the gain of the relay channel is reduced. Still, a gain of 6.3 and 9.5 dB can 
be obtained at FER = 10~ 4 for p = and 20 dB. 

The performance of C 4 i on the four-relay single-antenna channel is illustrated in Fig. [U The 
frame length is 16. For p = 20 dB, a gain of 21 dB (resp. 12.8 dB) at FER = 10 -4 is obtained 
with 4-QAM (resp. 64-QAM). 

5 A trivial suboptimal solution is to "turn on" the relay only when the p and SNR are high enough to give a better performance 
over the non-cooperative case. 
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B. Multi-Antenna Channel 

1) A (1, 2, 1)-Relay Channel: As an example for the case n r > n s , we consider the (1, 2, 1) 
channel. Here, we use the virtual iV-relay scheme with N = 2. As discussed in section Hll-A.21 
the diversity order of this scheme is comparable to the 2-relay single- antenna channel. To compare 
these two channels, we use the same code C 2 ,i- The performance is shown in Fig.|9j As compared 
to the 2-relay single-antenna channel, the (1, 2, 1) channel has a gain of 1.5 dB at FER = 1CT 4 
with p = OdB. With p = 20 dB, the two channels have essentially the same performance. In 
fact, the inter-relay cooperation in the virtual two-relay channel improves the receive SNR (3 dB) 
of the source-relay channel with antenna combining. Thus, the geometric gain p is increased 
effectively. However, as stated before, the global performace is not sensitive to p for large p's. 
This is why there is a gain only with small p's. 

2) A (2,2,2)-Relay Channel: In the case of n s > n r , we consider a single-relay channel 
with two antennas at each terminal. The code Ci j2 is actually the 4x4 Perfect code. For the 
non-cooperative scenario, we take the same code for fairness of comparison. More precisely, the 
non-cooperative channel we consider here is equivalent to a cooperative channel with 7r 3 = 
and 7Ti = 7r 2 . As shown in Fig.[T0l the gain of the cooperative channel over the non-cooperative 
channel is much less significant in the SNRs of interest. This is because the diversity order of 
the Perfect code in the two-antenna non-cooperative channel is already 4 and a diversity gain 
does not play an important role in the scope of interest. Note that at FER = 1CT 5 , the gain of 
the cooperative channel with p = dB over the non-cooperative channel is 2 dB for 4-QAM and 
3dB for 64-QAM. Also note that the difference between different p's is within ldB. 

3) A (2,2,2)-Relay Channel with Shadowing: In this scenario, we consider the shadowing 
effect of a wireless channel. Assume that each link between terminals is shadowed. Mathemat- 
ically, the channel matrix is multiplied by a random scalar variable, the shadowing coefficient. 
Suppose that this variable is log-normal distributed of variance 7dB [17] and that the shadowing 
is independent for different links. Fig. [TT] shows the performance of the cooperative channel 
with the use of Ci j2 (4-QAM) at p = OdB, as compared to the non-cooperative channel. FER 
is the averaged frame error rate on the channel fading and the shadowing. As shown in Fig.fTTl 
the slope of the FER-SNR curve of the non-cooperative channel is reduced as compared to 



February 1, 2008 



DRAFT 



22 



the non-shadowing case, in the scale of interest 6 . Since the shadowing is independent between 
different links, the cooperative channel mitigates the shadowing effect and we get a larger gain 
over the non-cooperative channel than in the non-shadowing case (Fig. |10(a)j )- At FER = 10~ 4 , 
this gain is 8dB, in contrast to 1.2 dB in the non-shadowing case (Fig. |10(a~ >. 

VIII. Conclusion and Future Work 

In this paper, a half-duplex MIMO amplify-and-forward cooperative diversity scheme is stud- 
ied. We derived the optimal diversity-multiplexing tradeoff of a MIMO Rayleigh product channel, 
from which we obtain a lower bound on the optimal diversity-multiplexing tradeoff of a MIMO 
NAF cooperative channel. Moreover, we established a lower and upper bound on the maximum 
diversity order of the proposed MIMO NAF channel and showed that they coincide when the 
numbers of antenna satisfy certain conditions. Based on the non-vanishing determinant criterion, 
we constructed a family of short space-time block codes that achieve the DMT of our MIMO 
NAF model. Our construction is systematic and applies to a system with arbitrary number of 
relays and arbitrary number of antennas. Numerical results on some explicit example codes 
revealed that significant gain in terms of SNR can be obtained even with some non-optimized 
parameters. This gain is much more important in the single- antenna case than in the MIMO 
case. Fortunately, in reality, it is also the case that we need cooperative diversity only when 
local antenna array is not available. 

Nevertheless, it still remains two important open problems to solve : 

1) optimization of the power allocation factors: Based on the statistical knowledge of the 
channel (notably p), how to choose the factors 7Tj's in order to optimize the code performance 
according to certain criteria? We set iti = 27r 2 = 2tt 3 for simplicity. It is clear that the 
optimal DMT is independent of these parameters. However, in practice, for different p and 
SNR, the factor 7r/s are significant for the performance {e.g., the error rate performance). 
How to analyze the impact theoretically is an interesting future work; 

2) optimization of the matrix B: we set B to be identity matrix and derived the lower 
bound (fT2l . Based on the receiver CSI at the relay, is there an optimal matrix B that 
gives a better DMT than the lower bound (fT2l ? This problem is independent of the code 

6 With shadowing, the slope converges very slowly to the diversity order of the fading channel. 
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we use. Solving this problem may lead to solution for the exact DMT of the MIMO NAF 
channel. 

Appendix I 
Preliminaries to the Proofs 

For sake of simplicity, we use the dot (in)equalites throughout the proofs to describe the 
behavior of different quantities in the high SNR regime. More precisely, 
• for probability related quantities, 

,. logpi logp 2 
Pi = P2 means hm - — = hm 



SNR^oologSNR SNR^oo logSNR' 

for mutual information related quantities, 

X\ = X 2 means lim - — = lim 



SNR^oologSNR SNR->oo logSNR' 

• for sets, 

«Si = S 2 means Prob {s £ S\} = Prob {s e S 2 } ■ 

> , < , 2 and C are similarly defined. 
Definition 3 (Exponential order [6]): For any nonnegative random variable x, the exponential 
order is defined as 

Z = - lim (24) 

S SNR^oologSNR 

We denote x = SNR~^. 

Lemma 4: Let X be a ^-distribution random variable with 2t degrees of freedom, the 
probability density function of its exponential order £ satisfies 

'SNR-°°, for£<0; 

.SNR - ^, for£>0. 

Let S be a certain set, £j's be independent random variables with ~ x| t , and P s = 
Prob j^, . . . ,£ N ) e <S j, then we have 



P s = SNR- d with d= inf 

(a-£jv)es+ f-' 



where 5+ = «Sf| RiV+ - 
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Appendix II 
Proof of Proposition Q] 

In the high SNR regime, the outage probability is [10] 

P out (r log SNR) = Prob {logdet (1 + SNRAA 1 ) < r log SNR} (25) 

with A = GH. Let us define C H = HW and W = AA. The entries of G and H being 
i.i.d. Rayleigh distributed, Ch and W\H are two central complex Wishart matrices [18]. Let 
pLi > ■ ■ ■ > jii> and Ai > • • ■ > A; > be the ordered eigenvalues of Ch and W, then we 
have [18], [19] 

ii / i 

/(/*) = G myl Y[ fx™' 1 Y[(n k - n P f exp I - ^ 

k<p \ fc=l 

A, 



fc=i 
i 



6XP 1 "7" 



fe=l k<p 

with ir^ n and G m>i being the normalization factors. Hence, the joint pdf of (A, fi) is 
/(A,/i) = /(M)/(A|/i) 

= Ci,m,n n ' (Afc — A p )(/ifc — /i p ) 

fc<p 



fc=l 



exp | - jMk I det 



fc=i 



Aj 

exp I 

fk 



where C/ ;m n is the normalization factor. Define a, = — log \/ log SNR and Pi = — log /Xj/ log SNR 
for i = 1, . . . , I. Then, we have 

i 

f(a,0) = Ci, m , n (logSNR) 2 ' JJSNR- (n - /+1)Qfc SNR- {m - n)/3fc 

k=l 

I 

■ Y[ (SNR- ak - SNR- Qp )(SNR- & - SNR"^) 



k<p 

■ exp 



exp (-SNR- ( ^~ ft) ) 



(26) 



^SNR- & | <ld 

k=l 

First, we only consider $ > 0, Vz, since otherwise, exp f— J^ fe SNR~^ fe ) decays exponentially 
with SNR [10]. Then, we can show that 

det [exp f-SNR- ( ^- ft)N )l = SNR-^=^<^-^ + , (27) 
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To see this, let us rewrite Dj = det 



exp ( -SNR- ( ^- ft) 



as 



A = e 



-E t SNR-C«i-*) det 



-SNR _ < a i~' 3 i'+SNR-( a !~' 3 i) 



-SNR"( a i-ft)+SNR-( a !-' 3 !) 



e -SNR-(-l-W det 



e -SNR-C-i-ft) det 



e -SNR~( Q i-' 3 i) _ g -SNR-( Q i-ft) 



g-SNR"^!"' 3 '-!' _ g-SNR-^i-ft) 
„-SNR-( Q i-' 3 !) 



e -SNR-("i-' 3 i) / j _ e -SNR-( a i-ft) 



g-SNR-^i-' 3 '-!) ^ _ g-SNR-^i-W 



J i j=l 

e -SNR~ (Q, i-i~' 3 i ) +SNR-( Q! !-' 3 i) 
e -SNR _(Q i-i-' 3 i ) +SNR-( Q! !-' 3 ;) 



g -SNR~ (Q! i-i~' 3l) _ -SNR-^t-i - ^' 3 



g-SNR-^i-i-^-i' _ g-SNR-^i-i"' 3 ;' 



-SNR - ^-! - ^ 
-SNR~ (Q i-i-' 3 i ) ^ _ 

-SNR-( a !-i-^-l) ^ 



,-SNR" (a ! 



-SNR~ (q 



l-l 



no 



— e 



-SNR-^i-^ 



i=i 



where the equations are obtained by iterating the identity SNR a ±SNR b = SNR a for a < 5. 
Since 1 - e~ x « x for x close to 0+, we have 1 - e -™R-(<*i-W = SNR~ (a * _A) if a, > pi 
and 1 — e~ SNR ( ° l 01 = SNR otherwise. As shown in the recursive relation above, we must 
have ctj > Vi, so that does not decay exponentially. In this case, we have e~ SNR <a * 01 = 
SNR°,Vi. Thus, we have D t = SNR~^ i < i( " l ~ ft)+ A-i, and in a recursive manner, we get d27b . 
Finally, we can write the outage probability as 



where 



P out (r log SNR) 



0(r)4 J( a> /3): - a k ) + < 



f(a,0)dadp 



0(r) 



r, „ . . „ , a* > Pi > 0, \n 



fc=i 



A < ■ ■ • < A 

is the outage region in terms of (a, /3) in the high SNR regime. Since SNR _afc — SNR~ ap (resp. 
SNR-^ - SNR~^) is dominated by SNR~ Qfc (myp. SNR~ A ) for k < p, from C3 and we 
have 

P out (r log SNR) = / SNfT'Wad/S 
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with 

i i i 

d a ,p = - i + + ~ n + 1 ~ + ~ ^) + - 

i=l i=l j=l i<j 

Let P out (r log SNR) = SNR -< ^ (r) . Then, we have 

d* A {r) = mfd a; p. (28) 

The optimization problem (l28t can be solved in two steps: 1) find optimal (3 by fixing ex, and 
then 2) optimize ex. 

Let us start from the feasible region 

< p x < Qfi < (3 2 < a 2 < ■ ■ ■ < A < a, 

z 

in which we have ^^^^("i — /%) + = 0. Note that the feasibility conditions require that /3j's 

can only move to their left in terms of their positions relative to the ctj's and that $ can never 
be on the left of (3j for i > j. Let bi denote the coefficients of in cf a)( g. The initial values 
of frj's are 6^ = m — n + l — i where > bf > ■ ■ ■ > bf . As long as bi is positive, $ 
should decrease (pass the <x, just left to it) to make the objective function d a ;i g smaller, with the 
feasibility conditions being respected. Each time a /?, passes a a.,- from right to left with j < i, 
bi decreases by 1. When bi = 0, $ should stop decreasing. Therefore, the optimal region is such 
that 

d = bf ] -b*= [min - 1, &< 0) }] + , 2 = 1,...,/. 

For 2 such that 6^ < 0, (3* = cti. For ft? ) > 0, (3* is «j_i_ Ci if i — 1 — q > 1 and otherwise. 
Note that &*'s are independent of ctj's and only depend on (m, n, Z). This is why we can separate 
the optimization problem into two steps. After replacing the optimal /3 in d a p and some basic 
manipulations, we obtain 



d, 



a 



sA ( , l + k + A\ , 
Y j iq+l-2k+ \a k + (9 + i + l-2A;) 

k=l ^ L - J / fc=Z-A+l 



ait 



= a J ex (29) 

where a k is non-negative and is non-increasing with k. Hence, the optimal solution is a k = 

l,k = s + 1, . . . , I and a k = 0,k = 1, . . . , s, from which we have d\(s) = ^2 S+1 a-k- For 
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s > I - A, 



d K s ) = J2 q + 1 + 1 ~ 2k 

k=s+l 

= (q-s)(l-s). 



(30) 



For s < I - A - 1, 



(9+l-2A:+ 1 + ^ + A ]o; fc + ^ (g + Z + 1 - 2fc) a fc 

fc— -Li ^ L - J / fc=Z-A+l 



fc=s+l 



/ 



(q-s)(l-s)- 



(l-A-s) 



(3D 



By combining (SDl and JITb. we get (fTTb . 



Appendix III 
Proof of Theorem [T] 

The main idea of the proof is to get lower bounds on the DMT by lower-bounding the 
mutual information of the channel defined by © and ®. Since the multiplicative constants 
have no effects on the DMT, for simplicity of demonstration, we will neglect them and rewrite 
S = I + PP t and 

r F 
TPH TF 



H 



The mutual information of the channel H is 

l(x; Hx + z) = log det (i + SNRHH^ . 
Lemma 5: Let y be the generalized inequality for matrices 7 , then 

and there exists a matrix B satisfying the power constraint © such that 

1 + A max (PP f ) = SNR°. 



(32) 



(33) 



' AyB means that A — B is positive semidefinite. 
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Proof: (l32t comes from the definition of £. (l33l can be shown by construction. Let us 



take 



VSNRB 4 y/c ■ min {A-^ (HJT) , 1} • I. (34) 

Then, the power constraint © is always satisfied with c = (SNR -1 + T^ip) 1 . Since c = SNR°, 
we have 

1 < 1 + A max (PP f ) < 1 + SNR _Qmax = SNR (35) 

where a max is the exponential order of A max (G t G) and is positive with probability 1 in the high 
SNR regime [6], [10]. ■ 
By lemma|5]and the concavity of logdet(-) on positive matrices, we have 

logdet (i + SNRM f ) > logdet (i + SHRHH^ > logdet (i + SNR(1 + A max (PP t )) _1 ™ t 



with 



H 



F 
PH F 



Therefore, with B in (EH), we have logdet (l + SMRHH^ = logdet (l + SHRHH^ . Assume 
that in the rest of the proof, we always consider B being in the form (l34b . Then we have 



J max = maxX(x; VSNRHx + z) > log det 1 1 + SHRHH 



(36) 



where B is the set of matrices B that satisfy the power constraint ©. Define M = 1 + SHRHH , 
we have 

r i + SNRPP f SMRFWP^ 
S N RPHF^ I + S N R (FF t + PHWP^ 



M 



det 



det (A) &et{D -CA-'B) 



Using the identity 

A B 
C D 

and some basic manipulations, we have 

det(M) = det (I + SNRPP f ) det (I + SNRPP f + SNRPHQH^) (37) 

where fl = I — SNRP f (I + SNRPF f ) -1 F is positive definite. By the matrix inversion lemma 
(I + LCR) 1 = I - L(RL + C-y'R, we have 

[1 = (I + SNRP+P)" 1 . 
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From (l36l) and d37t . we can obtain two lower bounds on X max . The first one is 

X max > 21ogdet(I + SNFLFF*), (38) 

whereas the second one is 

Xmax > logdet(I + SNRFF 1 ) + logdet(I + SNRPHQWP ] ) 
= logdet(I + SNRFtF) + logdet(I + SNRfiiTPtRH") 
= logdet(I + SNRF 1 ^ + SNRHTTH) 

> logdet(I + SNRfTPtPif). (39) 
Since in CSh min {A max (HW) , 1} = SNR°, from © and (03), we have 

X max > log det (I + SNRW&GH) . (40) 

The outage probability is 

, n , r m i • „ , f 21ogdet(I + SNRFF t ) < 2rlogSNR, 

< 2rl ° gSNR} £Pr0b logd!t(I + SNRirt?GffJ< 2H„gSNR 



= Prob 1 2 log det (I + S N RFF^) < 2rlogSNR 

■Prob|logdet(I + SURWGVH) < 2rlogSNR 
= SNR~ ( dF ^ +d oH( 2r )) 

where the second line follows from the independency between F and GH. 

Appendix IV 
Proof of Theorem [5] and Corollary [Tj 

A. Proof of Theorem [2] 

As in the case of the single-relay channel, we need two lower bounds on the mutual informa- 
tion. Since the mutual information of the A-relay channel is the sum of that of the A single-relay 
channels, these two lower bound can be obtained directly from d3*%t and d4*0l 



l(x; VSMRHx + z) > 2A log det (I + SNRFF 1 ") 

N 

I(x; VSNRHx + z) > log det (I + SNRfftGjG^ 



i=i 
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The outage probability is upper bounded by 



Prob | J(x; VSNRHx + z) < 2Nr log SNR 

f 2iVlogdet(I + SNR J F J F t ) < 2iVrlogSNR, 1 
~ r ° \ E 4 log det (I + S N RHppiHi ) < 2Nr log S N R J ' ( } 

Let us denote Zj = log det (I + SHRHlG\GiH i) and aij the set of exponential orders of the 
ordered eigenvalues of H\G\GiH \. The pdf of aij is p ai = SNR~ dai where from (l29b . is 
nondecreasing with respect to the component- wise inequality, i.e., 

d a >. > da. if a- b <*v (42) 



Let us define 



r N q 

O g (r) 4 {a,}f =1 : ££(l-a i)fe ) + <2iVr 

L i=i fc=i 

Oi(r) 4 Joi : ^(l-a iife ) + <rl. 



fe=i 



Then, the outage probability is 

JV 



Prob <j J^Z < 2iVrlogSNR J> = Prob {{a 4 }f =1 e 9 (r)} 

= SNR _drelay 



where rf re i a y is 



N 



drelay = inf V" d Qi 
/ N 

= inf 7 inf d a . 



,i=l 



inf ( Vd GiHl (2iV^r) ) 



with the second equality from the fact that the minimal elements lie always in the boundary 
when d42b is true. For GiHi identically distributed for all i, (fl~3t is obtained by the convexity 
of d GH . 
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B. Proof of Corollary \J\ 

For simplicity, we prove the particular case N = 1 here. For N > 1, same method applies. The 
lower bound is a direct consequence of theorem|2| The upper bound can be found by relaxing the 
half duplex constraint, i.e., H = T \PH F] with all matrices being similarly defined as before. 
Define H = [PH F}. First, since I y TT, we have 

X(x; VSNRHx + z) = log det (i + SNR#rT#) 

< log det (I + SNRFF t + SNRPHHT^ (43) 
= log det (I + SNRFF t + S N RG ( S N RBHWB^)G ] ) (44) 

< log det (I + S N RFF^ + S N RGG r ) (45) 

which means that the channel H is asymptotically worse than the channel [G F] in the high 
SNR regime. Thus, we have d F (0) + d G (0) > d NAF (0). 
Then, since I y T^P^PT, another bound is 

I{x; VSNRHx + z) = log det (I + SNRrFFT 1 + SNRTPHH f PT f ) (46) 

< log det (I + SNRTFFT 1 ) + log det (I + SNRTPHH^V) (47) 

< log det (I + SNR-FtF) + log det (I + SURH f H) (48) 

from which we have dp(0) + du(0) > d NAF (Q). 

When the channel is Rayleigh, proposition [l] applies. As indicated in remark 1, we have 
dQH = mm{dc, da} for \m — n\ > I — 1 and the lower bound and the upper bound in (fl4l 
match. 

Appendix V 
Proof of Theorem [3] 

To prove theorem[3l it is enough to show that in the high SNR regime, an error occurs with 
the rate-n NVD code X only when the channel is in outage for a rate -r. To this end, we will 
show that the error event set of X is actually included in the outage event set, in the high SNR 
regime. 
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A. Outage event 

For a channel H, the outage event at high SNR is [10] 

0(r) = {H : logdetfl + SNRM 1 ) < rlogSNR}. 

Let us develop the determinant as 8 

<? 

det(I + SNRM 1 ) = 1 + Sm l Di(HH f ) 

i=i 

where D^M) is the sum of (?) products of i different eigenvalues of M . In particular, we have 
D X (M) = Tr(M) and D„(M) = det(M). Let \ denote the 2 th smallest eigenvalue of HH and 
a,i denote the exponential order of \, i.e., \ = SNR - " 1 with a,\ > a 2 > • • • > a q . Then, we 
have 

A = SNR" E *=«-*+i a * fari = l,...,g 

since Y^l= q -i+i a k ^ s me smallest among all the combinations of % different a's. Now, we are 
ready to write 

0(r) = \h : 1 + SHRi D i( HH *) < SNRr 
I i=i 



{# : SNR* ACM 1 ) <SNR r , Vi=l,...,g} 



o; : « 



«fc ] < r, Vi = l,...,q 

\k=q—i+l 



1 



ol : a fc > (g - j) - r, Vj = 0, . . . , q - 1 

k=j+i 



(49) 



5. Error event of a rate-n NVD code 

Let us now consider the error event of a rate-n NVD code X. We will follow the footsteps 
of [11]. Using the sphere bound, the error event of ML decoding conditioned on a channel 

n 

8 To see this, consider the identity det(M — xl) = ( — 1)" Y\_( x ~ ^») wnere -V s ls tne eigenvalues of M. 
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realization H is 

dl ir 



E H C \ W : ||W||; > 



4 

= ju> : — w > 77 1 

where W is the AWGN matrix with i.i.d. entries; <i min is the minimum Euclidean distance 
between two received codewords, i.e., o^ in = min ||i7Ax || F ; w is the exponential order of 
||W||p (~ X2n R n T ) an d V i s mat °f V^min- Therefore, the error probability conditioned on H is 

Pe h <Prob{-w >r]}= SNR d£ « 

where by lemmaEl we have 

{inf n R riTW = 0, for 77 < 0; 
(50) 
00, for 77 > 0. 

Then the average error probability becomes 



^ = / Pe H PH(H)dH 
< [ PH (H)dH 

Jri<0 



>n<a 

= Probjr/ < 0}. 
Therefore, we get the error event in the high SNR regime 

St [a : r/<0}Cp|{a : f < 0} (51) 

with £ being any lower bound on 77. Using the same arguments as in [11], with a rate-71. NVD 
code, we can get q lower bounds on d^ in 

SNR" = < in (a) > SNR<^\ j = 0, . . . , q - 1 



with 



3V 71 ? + 1 ^ ? + 1 

i=q~j J 



Finally, from (BTT) and d52l . we get 

£ Q {<* ■ Sj < 0, Vj = 0, l} 



E - ~~ •?') — r ' v i = ° 5 1 

.k=j+l 
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which implies that 

dx(r) > d out QrJ . 

Appendix VI 
Proof of Theorem |4] 

Consider the channel A = diag(.Hi, . . . ,H N ) with Hi being similarly defined as H in © 
except that G,H,B in © are replaced by Gi,Hi,Bi, respectively. Since one channel use of A 
is equivalent to 2N channel uses of an A-relay NAF channel, i.e., 

Ca = 2AC NAFiJV (53) 

where Ca and C NAF>N are the capacities of the channel A and the equivalent A-relay NAF channel, 
measured by bits per channel use. Therefore, we have 

cWr) = <C,*(r) = dT(2Nr) (54) 

where the first equality comes from the fact that the outage upper bound of the tradeoff can be 
achieved [6] and the second comes from (l53l) and the definition of outage since 

P* t {2Nr) 4 Prob{C A < 2Nr logSNR} = Prob{C NAF>JV < rlogSNR} = P o N u T ,iV (0- 

On the other hand, by using a code C defined above, an equivalent channel model of the 
A-relay channel is 

Y = VSNRAX + Z 

with X E X. By theorem[3l we have 

d c (r) = d x (2r) > <(2Ar). (55) 

From (1541 and d531 ). we obtain 

d c (r) > d NAF , N (r). 
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Fig. 6. Two ways of extending Q(i) up to K. 



Appendix VII 

Cs IS NOT A NORM IN Q (( 8 , y/b) 

We prove, in this appendix, that Cs is not a norm of an element of K = Q (Cs, Vo). Assume 
that Cs is a norm in K, i.e., 

3x E K, N K/ Q i(j}) (x) = Cs- (56) 

Consider now the extensions described in Fig.|6j From d56l) . by considering the left extension 
of Fig.|6l we deduce that 

N K /Q(i)(x) = N Q((s)m) (N KM(s) (x)) =(s-r (Cs) = -%■ (57) 

Now, we deduce, from the right extension of figure |6] that 

N K /®(i){x) = N Q ^ m) (^k/q^VS)^)) = ~ l - (58) 
Denote y = N K ^^ ^\(x) E Q (i, V§). Then, the number z = 1+ 2 V ^ • y has an algebraic norm 
equal to i, and belongs to Q (z, v^)- In [12], it has been proved that i was not a norm in 
Q (i, \/5). So, Cs is n °t a norm in K. 

Appendix VIII 

Ci 6 IS NOT A NORM IN Q (Ci 6 , 
The proof is similar to the one of appendix IVII1 First, we assume that Ci6 is a norm in 

K=Q(Cie,>/5), i.e., 

3x e K, V K / Q ( Cl6 )(x) = Cie- (59) 
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We deduce that 

N K /Q(i)(x) = iV Q(Cl6)/Q(i) {N Kmw) (x)) = Cie • t (Cie) • r 2 (Cie) • ^ 3 (Cie) = -i- (60) 
But we also have, 

^K/Q(i)(x) = ^(i^/QCi) (^K/Q^VS)^)) = (61) 

Denote y = N K ^r ^(x) G Q (i, \/5). Then the number z = 1+ 2 V ^ • y has an algebraic norm 
equal to % and belongs to Q (i, \/5), which is a contradiction. 
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(a) 4-QAM 



Non-cooperative 
Coop, p = dB 
Coop, p = 5 dB 
Coop, p = 10dB 
Coop, p = 20dB 




24 36 

Receive E b /N (dB) 



, i i 

48 51 



(b) 64-QAM 

Fig. 7. Single-relay single-antenna NAF channel, Rayleigh fading, Golden code. 
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Receive £ b /No (dB) 



(a) 4-QAM 




12 24 36 48 54 



Receive E h /N (dB) 
(b) 64-QAM 

Fig. 8. Four-relay single-antenna NAF channel, Rayleigh fading, Caa- 
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Fig. 9. Single-relay (1,2, 1) NAF channel vs. two-relay single-antenna NAF channel, Rayleigh fading, Ci,\ with 4-QAM. 
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6 12 18 22 



Receive E b /N (dB) 



(a) 4-QAM 




(b) 64-QAM 



Fig. 10. Single-relay (2, 2, 2) NAF channel, Rayleigh fading, 4x4 Perfect code. 
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12 24 36 48 



Receive E h /N (dB) 



Fig. 11. Single-relay (2, 2, 2) NAF channel, Rayleigh fading, log-normal shadowing with variance 7dB, 4-QAM, 4x4 Perfect 
code. 
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